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PHYSICS OF STAR-PLANET MAGNETIC INTERACTIONS 


Antoine Strugarek 


Abstract. Magnetic interactions between a planet and its environment 
are known to lead to aurorae and shocks in the solar system. The large 
number of close-in exoplanets that have been discovered so far triggered 
a renewed interest in understanding magnetic interactions in other star- 
planet systems. Multiple magnetic effects were then unveiled, such as 
planet inflation or heating, planet migration, planetary material escape, 
and even some modifications of the host star apparent activity. Our 
goal here is to lay out the basic physical principles underlying star- 
planet magnetic interactions. We first briefly review the hot exoplanets' 
population as we know it. We then move to a general description of 
star-planet magnetic interactions, and finally focus on the fundamental 
concept of Alfvén wings and its implication for exosystems. 


1 Introduction 


1.1 The hot exoplanets population 


In the long trailing wake of the discovery of 51 Peg b (Mayor & Queloz 1995), 
more than 4000 exoplanets have been detected today. Space-based observatories 
such as Kepler and Corot have revolutionized this field by bringing a long-awaited 
wealth of exosystems, including unexpected massive planets on very short-period 
orbit (the so-called hot Jupiters). Complementing transit and radial velocity de- 
tection techniques, ground-based coronography further allowed the direct imaging 
of distant exoplanets and revealed the existence of planets on very wide orbits (for 
a review, see Pueyo 2018). The heterogeneous population of detected exoplanets 
is shown in Fig 1 as a function of the hosting star mass (vertical axis) and the 
orbital semi-major axis (horizontal axis), in units of stellar radius. The top panel 
presents the aggregated histogram of the known exoplanets as a function of their 
semi-major axis, which is heavily biased towards planets on compact orbits with 
a semi-major maxis typically smaller than 20 stellar radii’. Such close-in planets 


1Note that in the solar system, the semi-major axis of Mercury is about 83 solar radii 
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are expected to strongly interact with their host star through a variety of physical 
processes: they receive a strong irradiating stellar flux that affect their upper and 
lower atmospheres (Gronoff et al. 2020), they are subject to intense tides (Mathis 
2018), and they orbit in a strongly magnetized medium. We will focus here on the 
latter and refer the reader to the cited reviews for the other types of star-planet 
interaction. 


1.2 Observational status of star-planet magnetic interactions 


The search for star-planet magnetic in- 
teractions has been a continuous his- 4217 known exoplanets 
tory of positive and negative detec- ' 
tions. The first firm detection of a sig- 200 
nal that could be related to star-planet 
magnetic interactions was arguably re- 
ported by Shkolnik et al. (2003), fol- 


@ = Earth-like @ = Mercury-like 
T 


lowing a seminal theoretical argument e ee 
for the existence of star-planet mag- 9e * « 
netic interaction by Cuntz et al. (2000). € o 


They monitored at that time 5 promis- 
ing star-planet systems, and detected 
a modulation associated with the hot 
Jupiter orbital period in the K band of 
Ca II of HD 179949. The other sys- s 
tems did not show any direct evidence Semi-major axis (R.] 
of star-planet magnetic interactions at 

that time. Such signals have been ob- Fig. 1. Exosystem distribution from exo- 
served since for different systems (see planet.eu. The bottom panel shows the dis- 
the review by Shkolnik & Llama 2017) tribution as a function host-star mass (ver- 
such as famous compact exosystems  tical axis) and semi-major axis (horizontal 
like 55 Cne or HD 189733 (Cauley axis), colored by the planet mass (My sin i). 
et al. 2019). The star-planet interac- The top panel shows the same distribution 
tions tracers found in the stellar activ- only as a function of the semi-major axis 
ity indicators have always been subject normalized to the stellar radius. The clas- 
to a large variability with very clear gical habitable zone (Kopparapu et al. 2013) 
signals at some epochs, and no detec- is shown as a grey band for an Earth-like 
tion at others (see for instance Cauley planet. The P, = 10 days is labelled by 
et al. 2018; Shkolnik et al. 2008). We the magneta line. Finally, a ’magnetic influ- 
believe today that at least part of this ence limit’ line is shown in red (see 82.1). 
hide & seek stems for the intrinsic vari- 

able nature of star-planet magnetic interactions. Star-planet magnetic interactions 
are controlled by the large-scale magnetism of stars, that can be topologically com- 
plex (Donati & Landstreet 2009; See et al. 2015) and time-variable due to eruptive 
events (e.g. coronal mass ejection, on a monthly time-scale) or magnetic cycles 
(on the yearly to decadal time-scale). As a result, star-planet magnetic interac- 
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tion likely varies in nature and strength along planetary orbits as well as longer 
timescales. 

Historically, tracers of star-planet magnetic interaction have been searched for 
in radio. Indeed, compact exosystems look alike planet-satellite systems in our 
own neighbourhood. The interaction of Io in the magnetosphere of Jupiter is 
known to lead to cyclotron-maser radio emissions at the poles of Jupiter (Zarka 
2007). At the same location, strinkingly clear UV spots are present at the pole of 
Jupiter (Clarke et al. 2002a), tracing the footprints of the magnetic interaction of 
Jupiter with Ganymede, Io, and Europa. Until recently, no clear radio detection 
of star-planet magnetic interactions have been successfully obtained (for a recent 
review see Zarka 2018). In 2019-2020, two new observations hinted toward these 
long-awaited radio emissions. The first one was a cautious suspicion of a radio 
emission in the 7 Boo system by Turner et al. (2019). The second one was the 
detection of coherent radio emission from the red dwarf GJ 1151 by Vedantham 
et al. (2020). The latter is an indirect detection in the sense that the observed 
stellar radio emission is peculiarly intense compared to what one would expect from 
such a slowly-rotating star. The authors showed that the observed emission could 
be compatible with star-planet magnetic interaction with an Earth-like planet on 
a 1 to 5 days orbit, but we have not detected this planet so far. A firm detection of 
this hypothetical planet through transits or radial velocity would help strengthen 
our confidence in this detection scenario today. 


2 Magnetic interactions in compact exosystems: overview 


2.1 Sub- vs Super-alfvénic interaction 


Planets in close-in orbits are generally thought to have their orbit circularized by 
efficient tidal interactions (Souchay et al. 2013). We will therefore consider here 
only planets in circular orbit, for the sake of simplicity. The characterization of 
star-planet interactions requires to estimate the plasma conditions along the orbit, 
that are set by the stellar wind. We will not detail stellar wind theory here and 
refer the reader to the previous Evry Schatzmann school proceedings by Brun & 
Strugarek (2020)[Johnstone, this volume] where the basics of magnetized stellar 
winds are introduced (the interest read can start by the seminal papers of Parker 
1958; Weber & Davis 1967). Let us summarize it quickly as follows. Cool stars like 
the Sun possess an external turbulent convective envelope that is the seat of an 
efficient dynamo process. The associated large-scale magnetic field then sculpts 
the environment of the star. In addition, turbulent photospheric motions are 
thought to excite magnetic perturbations, likely in the form of Alfvén waves, that 
propagate along this magnetic field and ultimately deposit energy upper the stellar 
atmosphere. In the case of the Sun, this leads to a very hot corona reaching more 
than a million Kelvin. This extremely hot corona then leads to the existence of 
a stellar wind, accelerating outwards and reaching super-sonic and super-alfvénic 
speeds. At the Earth orbit, the solar wind speed averages between 400 and 700 
km/s, and reaches an alfvénic Mach number M, between 9 and 10! In the context 
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Fig. 2. Schematic of alfvénic Mach numbers in compact exosystems. The top panel 
illustrates the star (yellow), the planet (blue) and its orbit (black), and the direction 
associated with the alfvénic Mach numbers (Eq. 2.1-2.3). The bottom panels show the 
values of these numbers as a function of orbital distance. The left panel corresponds to 
a solar-like star, and the right panel to a young, fast-rotating Sun. 


of exosystems, one can then wonder what is the speed of the stellar wind around 
close-in planets. 

We can get a first estimate by computing a simple Weber-Davis (Weber & 
Davis 1967) wind solution under the influence of rotation. Several alfvénic Mach 
numbers can be computed, each characterizing different aspects of star-planet 
magnetic interactions. They are summarized in Fig. 2 and defined as follows 


Jp = y (2.1) 
UA 
; UU — Uke 
Masa = [ne — in| T wel, (2.2) 
miotal = |v" — Vicep| (2.3) 
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where v™ is the stellar wind velocity field, which we consider to be primarily 
in the radial star-planet direction (i.e. the accelerating component, denoted v?) 
and in the orbital direction (due to stellar rotation, denoted v7). The Kepler 
velocity Vkep is assumed here to be purely in the azimuthal direction (circular 
orbit). Finally, va = B, /./nuopw is the local Alfvén speed at the planet orbit based 
on the plasma properties in the stellar wind (where po is the vacuum magnetic 
permeability). 

In the bottom panels of Fig. 2 we illustrate the different alfvénic Mach numbers 
for a solar-like star on the left and a fast rotating Sun (e.g. a young Sun) on 
the right. The classical alfvénic Mach number, based on the radial velocity of 
the stellar wind (Eq. 2.1), is shown by the red dashed lines. For the Sun, this 
simplified 1D model predicts an average alfvénic point (M1242! = 1) around 17 
solar radii (i.e. 0.08 AU). The orbital Alfvénic Mach number is quite different as 
it is based on the relative orbital motion of the planet in the stellar environment. 
At first order, the corona of a star can be approximated to co-rotate with the star 
below the alfvénic point, and to rotate with a profile v7 œ r^? conserving angular 
momentum in the super-alfvénic region outside. 

The orbital alfvénic Mach number Mit?! is shown in blue in both panels. It 
is super-alfvénic only very close to the star, due to a very fast keplerian motion. 
We see that it remains nevertheless sub-aflvénic for most of the domain for a 
solar-like star. In the case of an extremely fast-rotating star (right panel), the 
co-rotation radius (vertical dashed line) can be in the sub-alfvénic domain of the 
stellar wind. At the co-rotation radius, M ^*^! = 0 by definition, which means 
that the magnetic interaction vanishes. Note that it is nevertheless a non-stable 
equilibrium of the star-planet system (see for instance Strugarek et al. 2014). 

One can conclude from this crude analysis that many of the planets orbiting 
below 20 R, are likely to be in a sub-alfvénic interaction regime, with all Alfvénic 
Mach numbers below 1. We illustrate this by applying our simple 1D Weber-Davis 
wind model to stars with varying mass and rotation rate (typically here between 
2 and 28 days). We obtain an estimate the Alfvénic transition, which appears as 
the red area along with the exoplanet distribution in Fig. 1. Note that some very 
close-in planets could still be in a super-alfvénic interaction state in their orbital 
direction, which still allows them to be magnetically connected to their host star 
in the radial direction. We now turn to the basic description of different magnetic 
connection scenarii in the sub-alfvénic regime. 


2.2 Different regimes of sub-alfvénic interactions 


A planet on a short-period orbit can be thought of as a perturber placed in a 
likely non-axisymmetric interplanetary medium. The perturbations it triggers can 
be decomposed on magneto-hydrodynamic (MHD) waves traveling away from the 
planet location (we will give the full description of this in 83.1). A net Poynting 
flux is therefore channeled away from the planet by waves, along what is often 
referred to as Alfvén wings (see § 3). These alfvénic perturbations can travel back 
and forth between the star and the planet due to reflection along the magnetic field 
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gradients. In the context of exosystems, a simple estimate of the ratio between 
the alfvénic travel time and the advective time-scale of the flow across the obstacle 
can be found in Strugarek (2018). In essence, what is found is that for obstacles 
constituted by a planet and its magnetosphere, alfvénic perturbations generally do 
not have time to perform a back and forth travel between the planet and the low 
atmosphere of the star during the advective crossing time-scale?. This situation 
corresponds to the so-called pure Alfvén wing case (Neubauer 1998). Other scenari 
can be achieved, depending on the propagation time of these waves and on the 
conductive properties of the planet. For instance, Laine & Lin (2012) showed that 
planets with no magnetosphere and high internal conductivities could lead to the 
opposite situation where Alfvénic perturbations of have time to perform multiple 
back-and-forth travels. The later case is generally referred to as the unipolar 
inductor interaction regime. A full description of the different interactions regimes 
can be found in (Strugarek 2018), we refer the interested reader to this review and 
the references therein. While Alfvén wings have been observed on Earth (e.g. 
Chané et al. 2012) and for natural moons of solar system's giants planets (e.g. 
Bonfond et al. 2017; Clarke et al. 2002b, and references therein), the distinction 
between the different scenarii has proved to be difficult observationally. In all cases, 
though, two Alfvén wings always form in the (Bẹ, v" — v**P) plane. Depending on 
the Alfvénic Mach numbers (82.1), both wings can connect onto the host star; only 
one may do so while the other extends away toward the interplanetary medium, 
or both may head away from the host star. As a result, the knowledge of the 
magnetic configuration between the host star and the orbital path of the planet is 
mandatory to assess star-planet magnetic interactions. 


3 Magnetic interactions in compact exosystems: the Alfvén wings 


We now turn to the theoretical concept behind our understanding of star-planet 
interactions: the concept of Alfvén wings (Neubauer 1980). 


3.1 Pre-requisites on magneto-hydrodynamic waves 


Let us recall here the properties of large-scale waves in magnetized plasmas. Ne- 
glecting rotation and stratification, three MHD modes can be identified: a slow 
MHD mode, a fast MHD mode, and pure Alfvén mode. The associated waves 
will be dubbed Alfvén wave (AW), slow MHD wave (SW), fast MHD wave (FW). 
Note that the term magnetosonic is generally used for purely perpendicular waves, 
we thus prefer to keep here the generic MHD denomination for these modes and 
waves. The derivation of the dispersion relation of these modes can be found in 
most good MHD textbooks (e.g. Goedbloed et al. 2019 to cite only one). We give 
only here the result and let the reader dig into such textbooks for their derivation. 


?Note that in general, the reconnection timescale is short enough to not play major role here 
in setting the type of star-planet magnetic interaction. It can nevertheless play a role in setting 
the absolute amplitude of these interactions. 
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We will write k the propagation vector of the wave, and w its frequency. We 
will assume that at the studied scales the magnetic field can be considered to 
be homogeneous and we write Bo = Bob. The angle between k and Bo will be 
written 0. Finally, the classical Alfvén and sound speeds are defined as 


B 
VA = 2 d (3.1) 
v HOPO 
^J Po 
é Ex ait. 3.2 
Po ( ) 


where pp the plasma density, Pj the plasma pressure and y the ratio of specific 
heats. The dispersion equation for a pure Alfvén wave is 


n = vAcosÓ, (3.3) 


and the dispersions for fast (--) and slow (-) MHD waves are 


A 
w| E 
NY 
N 

Il 


1 
2 | (e 472) x JJ (cà +02)? — Ac2v2, cos? 6 
2 s A s A s"A 


(8-1) + "IE + D — 4 cog? ^ (34) 


where B is related to the plasma ( through 


A 
2 


d 
7 _ y2uofo _ y 
p= = lS L lg. (3.5) 


A graphical representation of Alfvén waves can help guide our intuition on 
their properties. In particular, Friedrichs diagrams can be built for these three 
fundamental waves. For a given frequency, we represent in Figure 3 the phase speed 
of the wave vp = (w/ k)k for all propagation directions k, where k is the unit vector 
aligned with k. In each panel, the magnetic field direction b is represented by the 
vertical direction. The perpendicular direction e, is represented by the abscissa. 
In 3D, the Friedrichs diagrams are invariant by rotation along b. Asa result, we 
represent these diagrams in 2D, bearing in mind that the full characteristics can 
be obtained by rotation along b. 

We represent the phase speeds for the AW (black), the SW (blue) and FW 
(red). The AW and SW are degenerate for k — 0 and fall back on what is some- 
times called an entropy wave. Two illustrative cases are shown in Figure 3. On 
the left panel, 8 < 1 and the AW and FW have the same phase speed when k || b. 
On the right panel, B > 1 and the AW and SW have the same phase speed when 
k || b. We note that as @ increases, the sounds speed overcomes the Alfvén speed 
and the AW and SW converge essentially to the same the phase speed portrait as 
the dispersion relation of SWs converges to the one of AWs. 
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Friedrichs diagram: phase velocity 


Fig. 3. Friedrichs diagrams of the phase velocity of Alfvén waves (AW, black), fast MHD 
waves (FW, red) and slow MHD waves (SW, blue). The vertical axis is aligned with an 
ambient magnetic field, the horizontal axis represents any perpendicular axis. On the 
left panel the Alfvén speed is larger than the sound speed, on the right panel it is the 
opposite. 


The same diagram can be build for the group velocity of the waves. The group 
velocity is defined as 


(3.6) 


The first term corresponds to the parallel component of the group velocity (which 
is equal to the phase speed), and the second to perpendicular component of the 
group velocity. In the case of AWs, we directly obtain that only the parallel 
component remains and obtain 


vg =vab. (3.7) 


We immediately note an important result: the group velocity of pure AWs is 
degenerate as it does not depend on the wave vector k. As a result any wave packet 
of pure AWs has to propagate along the magnetic field b. The same calculation 
can be followed for SWs and FWs, and one obtains 


^ o cos Ô sin 0 
vt =v k âl, 3.8 
d P 2V1 — ø cos? 8 [1 + v1 — e cos? 0] (8) 
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Friedrichs diagram: group velocity 


Fig. 4. Friedrichs diagram for the group velocity of large-scale MHD waves. The layout 
is the same as in Fig. 3. Note that the diagram for Alfvén waves consists in the two 
black dots due to the degeneracy of their group velocity. 


PS 2 
where we have introduced o = 40/ (8 + 1) , and ni is the direct unit vector 


perpendicular to k in the (b, k) plane. 

'The group velocity indicates the propagation direction and speed of wave pack- 
ets. It thus constitutes an important information to characterize how energy is 
channelled by wave ensembles. They are represented for AWs, SWs and FWs in 
Figure 4 with the same layout as in Figure 3. FWs are completely unfocused 
as their group velocity can take any direction. SWs have a more focused group 
velocity around b. This focusing is nevertheless less and less tight as 8 increases. 
Finally, AWs are completely focused (black dots in Figure 4), as we already noted 
earlier. This will have major consequences on the existence of Alfvén wings, which 
we now turn to. 


3.2 The concept of Alfvén wings 


We will now apply the basic MHD concepts of Section 3.1 to the case of star- 
planet magnetic interactions. We want to characterize the stationary waves that 
can exist due to the existence of an obstacle in a magnetized flow. In this simplified 
view, the obstacle can be a conducting planet, or the ensemble of a planet with 
its magnetosphere. The detailed shape of the obstacle itself does not matter for 
the moment. 

We will consider the frame at rest with the obstacle. Let us define vo the flow 
velocity in this frame, such that the obstacle is subject to a stream v, — voh that 
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Fig. 5. Friedrichs diagrams of the phase velocity involving a Doppler-shift associated 


B = 0.04 
MA = 0.80 


with a large-scale flow impacting the obstacle. The Doppler-shift is labelled by the red 
circle. The direction of the incoming flow h is given by the red line. The two panels 
show the same interaction case with a different impact angle O of the flow. 


deviates from being perpendicular to b by an angle O. In this frame, a stationary 
wave solution requires that 
w=k-vo, (3.9) 


that is naturally the Doppler shift in this frame. If on applies the same reasoning 
for Doppler shift that the one for Alfvén waves developed in previous section, 
we remark that such waves are described in the Friedrichs diagram by a circle of 
radius v9 /2 centered on (— cos © v9/2, — sin © vp /2) in the (b, e, ) plane (again, in 
3D this is actually a sphere by rotation along h). We can intersect this circle with 
the three Alfvén waves characteristics in the phase velocity Friedrichs diagram 
as shown in Figure 5. The possible stationary waves then depend on the most 
important control parameter for star-planet magnetic interactions: the (total) 
Alfvénic Mach number, defined here as 


dM. (3.10) 
UA 

On the left panel of Figure 5, we consider the canonical case O — 0 for a 
sub-alfvénic interaction (MA = 0.8). On the right panel, we illustrate the effect 
of the flow inclination angle O. The Doppler-shift is labelled by the red circle. 
We rightfully note that it actually intersect only the characteristics of the SWs 
and AWs. If the interaction was super-alfvénic (MA > 1), the red circle would 
actually intersect the three characteristics. The Doppler, SW and AW circles all 
cross in the center. We find again the 'entropy wave', for which the wave vector k 
is perpendicular to vo and the group velocity is equal to vo. These waves form a 

queue behind the obstacle, in the direction of the impacting flow. 
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Fig. 6. Schematic of the concept of Alfvén wings. The obstacle, illustrated by the 
blackened areas, can consist of a planet and its hypothetical magnetosphere. A front 
(left panel) and side (right panel) are shown. The two wings are labelled in red and blue. 


We have seen in Section 3.1 that the group velocity of AWs is degenerate. 
As a result, thanks to the Friedrichs diagram analysis we know that pure AW 
wave packets will propagate along the so-called Alfvén characteristics (blue lines 
in Figure 5) defined as 


Cy = Vo EVA. (3.11) 


Finally, the case of SWs is interesting. Alike the pure AWs, they show an 
intersection with the Doppler shifted red circle in Figure 5. Nevertheless, their 
group velocity is non-degenerate. As a result, any excited SW wave packets will 
occupy a larger and larger volume as we move away from the obstacle. The energy 
and momentum transport by these waves will therefore be unfocused and are 
generally ignored in the context of star-planet or planet-satellite interactions. 

We can sketch the characteristics of sub-alfvénic star-planet interaction de- 
duced from this simple Friedrichs diagram analysis, as shown in Figure 6. The 
Alfvén wings always live in the (b, h) plane. The inclination angle between the 
Alfvén wings and the ambient magnetic field, O4, can be geometrically derived 
from Figure 5 to obtain 


Ma, cos OQ 


sinOy = à 
^' Vic M3 — 2M, sin 


(3.12) 


To summarize, a first analytical description of Alfvén wings yields the following 
two robust conclusions: 
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Fig. 7. Alfvén wings in the Kepler-78 system (Strugarek et al. 2019). The left panel is 
an illustrative view of the corona of Kepler-78 deduced from a 3D, MHD modelling using 
the observed magnetic map of Kepler-78 (Moutou et al. 2016). The right panels show the 
projection of the two Alfvén wings on the ecliptic plane for four different orbital phases 
orb. Figures are adapted from Strugarek et al. (2019) with permission of the AAS. 


e Only pure Alfvén waves and slow MHD waves are excited by the obstacle 
in the subsonic flow. Between the two, only pure Alfvén waves packets 
propagate in a focused manner, which leads us to consider only this wave 
population in the energetic discussion hereafter (§3.3) 


e The overall shape of the Alfvén wings does not depend on the specifics 
of the obstacle. It is fully determined by the flow impacting the obstacle 
and the ambient magnetic field in the flow. 


The latter has a particularly interesting consequence: one can predict the shape 
of Alfvén wings solely based on the structure of a stellar corona and the orbital 
motion of an exoplanet. We show the application of the Alfvén wings theory to 
the Kepler-78 exosystem in Figure 7 (for more details see Strugarek et al. 2019). 
On the left, a polytropic wind model is used to derive the plasma properties of 
the corona of Kepler-78, based on observed Zeeman-Doppler-Imaging magnetic 
map observed by Moutou et al. (2016). The orbit of Kepler-78b is shown by the 
dashed white circle. On the right, we apply the Alfvén wings theory to predict 
how wings change along the orbit of the planet. We show the projected trace of 
north (in red) and south (in blue) wings on the ecliptic plane. We observe that as 
the planet orbit the AW shape and sub-stellar points (red and blue star symbols) 
dramatically change due to changes in magnetic connectivity in the stellar corona. 
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3.3 Energetics of Alfvén wings 


Now that we have understood how to characterize the shape of Alfvén wings, let 
us turn to their amplitude. How much energy is involved in magnetized star- 
planet interactions? The energy flowing through the wings is characterized by the 
Poynting flux S defined as 


ExB 
=, (3.13) 
Ho 
Assuming ideal magneto-hydrodynamics, the electric field is simply E = —v x 


B. The projection of S along the Alfvén characteristics in the plasma reference 
frame is therfore 


T B2 
Saw =S- A. 2 anges. (3.14) 


[5 Ho 


Equation 3.14 gives a robust, maximum estimate of the power involved in the 
magnetic interaction. It corresponds to the case where the perturbations propa- 
gating in the wings are sufficiently large so that the incoming flow is completely 
halted inside the wings. 

In general, the exact Poynting flux depends on two important aspects of star- 
planet magnetic interaction: (i) the details of the interaction inside the grey area 
in Fig. 6: what are the plasma, magnetic, and conductive properties of the obsta- 
cle? and (ii) the size of the effective obstacle. The latter determines the thickness 
of the Alfvén wings, over which one has to sum to estimate the total power in- 
volved in the interaction. We illustrate these aspects in Fig. 8 where we show 3D 
simulations for three different magnetic topologies (Strugarek et al. 2015). The 
Alfvén wings are highlighted by the volume rendering of the currents in blue/red 
(negative/positive). All three cases are equivalent: in the two first cases only the 
orientation of the planetary field is reversed, while in the third case the stellar 
magnetic field is changed from a dipole to a quadrupole of similar amplitude. We 
strikingly see in these renderings that the relative topology of the magnetic field 
has an important influence on the strength of the wings, while their overall shape 
is control by the stellar atmosphere. 

Saur et al. (2013) have developed an analytical model to calculate the total 
Poynting flux over Alfvén wings, assuming the obstacle is magnetized. They found, 
for simple geometries, 


P= | S-dA x~ Ra? (1+ M2 - 2MAsinO) " Saw, (3.15) 
AW 


where R is the effective radius of the obstacle, and & quantifies the conductive 
properties of the obstacle?. These results were later confirmed and extended to 


35 = Xp/(Xp +a), where Xp is the Pedersen conductance in the planets ionosphere and 
Da = (1/uovA)(1 + M2 — 2M4 sin 9)-!/? is the Alfvén conductance. 
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Fig. 8. Illustration of the effect of the magnetic topology on Alfvén wings, from the 
3D numerical simulations of Strugarek et al. (2015). The two left cases are the same 
except for the orientation of the planetary field (grey arrow), showing a modulation of 
the wings in amplitude (the red/blue volume renderings show the current system of the 
Alfvén wings). The right panel is the same except for the stellar magnetic field that is 
quadrupolar, showing a change in the shape of Alfvén wings. Figures are adapted from 
Strugarek et al. (2015) with permission of the AAS. 


more complex geometries through 3D numerical simulation in Strugarek (2016, 
2017). The maximal power is obtained for high conductance in planet's atmop- 
sheres for which à = 1. In the solar system, à varies from a few 10^? to about 
1 for the different moons of Jupiter and Saturn (Kivelson et al. 2004; Saur et al. 
2013). The conductive properties of exoplanets are nevertheless largely unknown 
as of today. They depend for instance on the complex photo-chemical reactions in 
their upper atmosphere, and hence on both the irradiation recevied by the planet, 
and the planet atmospheric composition. Both are expected to change significantly 
along stellar evolution. Ab-initio calculations are today being performed to charac- 
terize (among other aspects) the conductive properties of exoplanets' atmosphere 
(e.g. see the Kompot code, Johnstone et al. 2019). 

Assuming a maximized interaction, the total Poynting flux P requires the 
knowledge of the plasma environment around a planet along its orbit. Saur et al. 
(2013) considered a simple Parker-like wind and applied their formalism to the en- 
semble of known exoplanets at that time. The largest uncertainty lies here in the 
lack of knowledge of the key stellar properties (rotation period, magnetic field) in 
the sample. For this review I revisited the updated ensemble of Kepler-observed 
exosystems while using the wind-modelling strategy developed by Ahuir et al. 
(2020). This strategy can be applied for the stars with a detected rotational pe- 
riod (McQuillan et al. 2014): a stellar Rossby number can then be estimated based 
on the scaling law derived from 3D numerical simulation by Brun et al. (2017). 
The stellar magnetic field and the subsequent wind is then computed following 
the approach proposed by Ahuir et al. (2020) that satisfies all the observational 
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contraints on stellar wind available to date. The resulting estimated Poynting flux 
(Eq. 3.15) assuming a planetary surface field of 4G for these systems is shown in 
Fig. 9. In spite of the large scatter due to the large variety of star-planet systems 
considered here, we see a clear overall trend P œ p. ? for the compact exosys- 
tems in the Kepler field. We highlight two populations in Fig. 9: exosystems in a 
very likely sub-aflvénic interaction regime (red circles), and exoplanets in a likely 
sub-alfvénic interaction regime (green squares). The former population is found 
to be in a sub-alfvénic interaction within all the error-bars of the models and ob- 
servations used here (see Ahuir et al. 2020 for more details on these error-bars). 
The latter population (green squares) is predicted to be sub-alfvénic only for part 
of our estimated magnetic field and coronal properties of the central star. It is 
nonetheless instructive to observe that sub-alfvénic interactions are highly unlikely 
to occur for planets with an orbital period longer than five days. 

The Kepler-70 system stands out in this diagram. It is composed of a hot 
sub-dwarf around which a super-Earth (Kepler-70b) orbit with an extremely short 
period. Our estimates predict a moderate stellar magnetic field, which translates 
into a quite weak star-planet magnetic interaction. As a result, small orbital 
periods do not necessarily imply strong star-planet magnetic interactions. This 
is also illustrated by the large vertical spread in Fig. 9: planets with the same 
orbital period can lead to absolute Poynting fluxes varying by more than 3 orders 
of magnitude (e.g. for P44, ~ 2 days) depending on the spectral, rotational, and 
magnetic properties of their host star. 

The case of Kepler-78 is also highlighted in the top left of Figure 9. The 
Poynting flux estimate was carried out using our empirical scaling law with no a 
priori knowledge of the stellar magnetic field. It nonetheless provides an estimate 
that is very close to the one obtained when carrying out fully 3D MHD simulations 
of Kepler-78, based on an observed magnetic maps with Zeeman-Doppler Imaging 
(see Fig. 7 and Strugarek et al. 2019). 

We have illustrated the theoretical derivation of the energetics of star-planet 
magnetic interaction with the Kepler sub-sample of exoplanets. This sub-sample 
is particularly interesting because asterosismology allows us to have access with 
good precision to the fundamental stellar parameters required to properly assess 
magnetic interactions. Many other exosystems as for instance HD 1897333 or 55 
Cnc (Cauley et al. 2019) also present promising configurations for sub-alfvénic in- 
teractions. Dedicated, 3D models of such systems are today required to properly 
assess both the energetics and the ephemerids of star-planet interactions (Stru- 
garek et al. 2019). One can hope that the coming years will see the definitive 
confirmation of star-planet magnetic interactions detection thanks to a combina- 
tion of advanced 3D models supporting this interpretation of peculiar observational 
features correlated with exoplanets' orbital periods. 


3.4 Angular momentum transfer and the secular evolution of stellar systems 


The mere existence of Alfvén wings implies an energy transfer that also contains 
in part an angular momentum transfer. Indeed, an orbiting planet interacts with 
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Fig. 9. Estimated Poynting flux power (Eq. 3.15) for the Kepler sub-sample of McQuillan 
et al. (2013) as a function of the orbital period of the planet. The red dots label the 
planets that are always predicted to exhibit star-planet magnetic interactions within all 
the error-bars of the model. The green squares label the planets that are predicted to 
exhibit such interactions only for the most optimistic error-bars. 


its environment and therefore is subject to an equivalent drag force. For the solar- 
system planets, the stellar wind is not quite dense enough for this force to matter 
at all in their orbital evolution. But it may not be systematically the case for 
close-in exoplanets orbiting in a dense and strongly magnetized stellar coronae. 


The magnetic torque felt by close-in exoplanets was systematically parame- 
terized through numerical simulations in Strugarek (2016, 2017). We show the 
example of one such simulation in the left panel of Fig. 10, for one particular 
magnetic topology (see Ip et al. 2004; Strugarek et al. 2015 for more details on the 
effect of topolgy). A clear magnetic torque originating from the magnetic tension 
around the planet (green line) produces most of the total net torque (grey line) 
felt by the orbiting planet. The main unexpected finding of these studies was that 
this net magnetic torque could make close-in planets migrate on a time-scale of 
hundreds of millions years. Henceforth, magnetic torques could in principle play 
a role in shaping the population of close-in exoplanets. This situation typically 
arises mainly for strongly magnetized dwarf stars, or early in the evolution of a 
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solar-like star. For main-sequence slow-rotating stars, magnetic torques play a 
more negligible role: for instance, the fastest migration timescale found in the 
Kepler sub-sample studied in Fig. 9 is predicted for Kepler-828 with a timescale 
of approximately 2 Gyrs. 


Nevertheless, it is important to consider magnetic torques for ensemble evolu- 
tion of compact exosystems, particularly during their early evolution. A system- 
atic comparison between tidal and magnetic torques can be carried out thanks to 
the developments in both theoretical fields over the last decade. We have carried 
out such a comparison in Strugarek et al. (2017b) and showed that in general 
both physical phenomenons add up to make orbital migration more efficient. We 
found that magnetic torques were likely dominating the migration process for 
fully-convective stars and young stars, while tidal torques were likely dominating 
otherwise. We summarize this instantaneous comparison in the right panel of Fig. 
10, for a typical TTauri star - hot Jupiter system (alike Tap-26, see Yu et al. 2017). 
We see that the torque can be dominated by either tidal torques (red regions in 
the top panel) or magnetic torques (blue regions) in a Pyot-Porn diagram. The 
bottom panel shows the migration timescale resulting from the sum of tidal and 
magnetic torques (see Strugarek et al. 2017b for more details). 


It is today mandatory to go beyond the instantaneous approach to understand 
how compact star-planet systems jointly evolve. Any orbital migration corresponds 
to an angular momentum exchange with the star, that ultimately affects the stellar 
rotation period, and thus its dynamo action and the sub-sequent evolution of the 
system. This is a challenging endeavour as it involves the complex stellar secular 
evolution along with non-linear physical phenomena such as stellar winds, and 
tidal and magnetic star-planet interactions. A seminal exploration was carried 
out simultaneously by Bolmont et al. (2012); Penev et al. (2012); Zhang & Penev 
(2014) exploring the angular momentum budget of a two-body star-planet system 
along stellar evolution, including tidal interactions. Recently, Gallet et al. (2018) 
followed a similar approach based on state-of-the-art stellar evolution tracks pro- 
duced by the STAREVOL code (Amard et al. 2016), and advanced stellar wind 
(Matt et al. 2012) and tidal interactions (Bolmont & Mathis 2016; Ogilvie 2013) 
modelling. In this study they focused on the hot exoplanets population predicted 
by such models for young open clusters such as Pleiades. A similar approach 
was also used with the ESPEM code (Benbakoura et al. 2019) to predict which 
star-planet systems were the most likely to undergo a coupled orbital-rotational 
evolution phase. Star-planet magnetic torques yet have to be included in these 
theoretical explorations. This tool can ultimately be used to produce synthetic 
exoplanet populations, as we recently demonstrated in Ahuir et al. (2021). 


Such studies are today highly valuable in the context of Kepler and TESS (and 
in the future PLATO) as the number of the observed exosystems keeps growing 
fast. They will allow to disentangle the relative importance of the different physical 
effects in shaping the population of exoplanets thanks to dedicated, careful cross- 
comparisons between observed and synthetic samples. 
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Fig. 10. Illustration of magnetic torques due to star-planet magnetic interactions. The 
left panel shows an angular momentum budget around the orbiting exoplanet from 3D 
numerical simulations (Strugarek et al. 2015). The main magnetic torque is dominated 
by the magnetic tension contribution in green. The right panel shows the comparison 
between tidal and magnetic migration torques in a (Prot, Porb) diagram for TT-Tauri — 
hot Jupiter exosystem (see Strugarek et al. 2017b for more details). The bottom right 
panel shows the total (tidal+magnetic) migration timescale in a logarithmic color-scale 
for the same system. Figures are adapted from Strugarek et al. (2015) and Strugarek 
et al. 2017a with permission of the AAS. 


4 Conclusions 


We have presented here the basics behind our physical understanding of star-planet 
magnetic interactions. Several grey areas nevertheless remain to be explored and 
understood. Among them, I would like to highlight in particular that we have 
almost not discussed the details of the obstacle itself, i.e. the planet and its 
hypothetical magnetosphere/ionosphere/etc... While the toplogy of star-planet 
magnetic interactions is largely controlled by the star itself, their strength is sen- 
sitive to the detailed properties of the obstacle. What are the conductivites in the 
different layers of the atmosphere and interior of exoplanets? This question cannot 
be trivially anserwed today, as exoplanets can have wildly different chemical reac- 
tions compared to solar system planets, thay can receive a much different ionizing 
flux from their star, and they can be subject to much stronger tidal interactions 
with their host star and/or other planets in the system. 

In spite of these hurdles that have to be overcome on the theoretical side, the 
characterization of star-planet magnetic interactions should allow us to probe the 
magnetism of exoplanets. This grand goal is worth such efforts, as today we can 
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only characterize (with not so much precision) the magnetism of the few planetary 
bodies in the solar system. One may await a revolution in our understanding of 
planetary magnetism if the hunt for star-planet magnetic interactions succeed in 
the future. I hope this review will help motivate more and more young researchers 
to dive into the fantastic open world of star-planet interactions. 


I warmly thank J. Saur and K. Khurana for useful discussions on planet-satellite and star-planet 
interactions, and the exoplanet Team behind the exoplanet.eu website. I want to thank also J. 
Ahuir for a careful read of this manuscript. I am grateful to INSU/PNP for its financial support 
for this research. The numerical simulations shown in this manuscript were made possible thanks 
to GENCI time allocations over the years. 
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